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(5) Verfahren zur Klassifikation einer meSbaren Zeitreihe, die eine vorgebbare Anzahl von Abtastwerten 

aufweist. beispielsweise eines elektrischen Signals, durch einen Rechner und Verwendung des Verfahrens 



) Bai dem Verfahren warden fur die Abtastwerte bedingte 
Entropion ermitteit, anhand derer ein InformationsfluS fur 
eine vorgebbare Anzahl zukGnftiger Abtastzeitpunkte be- 
stimmt wird. Anhand des Informationsflusses wird eine 
Klassifikation der Zeitreihe durchgefuhrt. Der Informations- 
fluS spiegelt nichtlineare Korrelationen zwischen den Ab- 
tastwerten wieder. Somit ist eine Klassifikation zwischen 
Zeitreihen, deren Abtastwerte nichtlinear korreliert sind und 
Zeitreihen, deren Abtastwerte stochastisch unabhangig sind, 
mdglich. 
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Beschreibung 

Die Erfindung betrifft ein Verfahren zur Klassifikation einer meBbaren Zeitreihe, die eine vorgebbare Anzahl 
von Abtastwerten aufweist, insbesondere eines elektrischen Signals, durch einen Rechner. 

Die Erfindung findet Anwendung in technischen Gebieten, in denen es von Interesse ist, aus geraessenen 
Zeitreihen auf das zukQnftige Verhalten der Zeitreihen zu schlieBea Diese Vorhersage des zukfinftigen "Verhal- 
tens* der Zeitreihe erfolgt unter der Annahme, daB die Zeitreihe nichtlineare Korrelationen zwischen den 
Abtastwerten der Zeitreihe aufweist. . . • t 

Erhebliche Bedeutung erlangt dieses Problem auch in verschiedenen Gebieten der Medizin, beispielsweise in 
der Kardiologie. Speziell in dem Problembereich des Pldtzlichen Herztodes kann es lebenswichtig sein, Fruh- 
warnzeichen des Pldtzlichen Herztodes zu erkennen, urn so frtth wie mdglich GegenmaBnahmen gegen das 
Eintreten des Plotzlichen Herztodes einzuleitea 

Es ist bekannt, daB eine Zeitreihe eines Elektrokardiogramms, welches nicht korreliert ist, ein nicht gef ahrde- 
tes Herz bezQglich des pldtzlichen Herztodes beschreibt Ein gefahrdetes Herz bezQglich des pldtzlichen 
Herztodes wird durch eine Zeitreihe des Elektrokardiogramms beschrieben, welches nichtlineare Korrelationen 
zwischen den Abtastwerten der Zeitreihe aufweist [1} Weiterhin ist es aus [1] bekannt, aus der graphischen 
Phasenraumdarstellung (Fourier-Transformation) zweier aufeinander folgender Herzschl&ge Zeitreihen eines 
Elektrokardiogramms zu ermitteln, die Herzen beschreiben, die bezQglich des Pldtzlichen Herztodes gefihrdet 
sind. 

Das in [1] beschriebene Verfahren weist alle Nachteile auf, die empirische Verfahren in sich bergen. Hierbei 
sind insbesondere die Fehleranfalligkeit graphischer Deutungen durch einen Menschen, das Problem des Set- 
zens einer Schranke, ab der eine Zeitreihe als gefthrdet klassifiziert wird, sowie Ungenauigkeiten in der 
Darstellung der Fourier-Transformierten auf dem Bildschirm als Nachteil des bekannten Verfahrens zu betrach- 
ten. 

25 Weiterhin sind Verfahren zur Bestimmung stochastischer, bedingter Entropien bekannt [2 J [3J #> 

Aus [4] ist ein Verfahren bekannt mit dem der zeitliche Verlauf der lokalen Sauerstoffspannung des Gehirns 
(tip02) ermitteit werden kann. 

Aus dem Dokument [5] ist ein Verfahren und eine Anordnung zum Vergieichen von Wellenformen von 
zeitveranderlichen Signalen bekannt 
30 Der Erfindung liegt das Problem zugrunde, ein Verfahren zu schaffen, urn eine Zeitreihe, die eine vorgebbare 
Anzahl von Abtastwerten aufweist, mit Hilfe eines Rechners schneli und verlaBlich zu klassifizieren. 
Das Problem wird durch das Verfahren gemaB Patentanspruch 1 geldst 

Bei dem erfindungsgemaBen Verfahren werden fQr eine vorgebbare Anzahl von Abtastwerten bedingte 
Entropien ermitteit Aus den bedingten Entropien wird ein InformationsfluB fiir eine vorgebbare Anzahl zukunf- 
35 tiger Abtastpunkte bestimmt, anhand dessen die Zeitreihe klassifiziert wird 

Weiterbildungen des erfindungsgemaBen Verfahrens sowie dessen Anwendungsmdglichkeiten ergeben sich 
aus den abhangigen Ansprfichen. 

Durch das Verfahren gemaB Patentanspruch 5 ist es mdglich, die KJassifikation zu beschleunigen, da nur 
anhand der Form des Graphen des Informationsflusses eine binare Klassifikation durchgefuhrt werden muB. Die 
40 Unterscheidung der Zeitreihe in einen ersten Zeitreihentyp und in einen zweiten Zeitreihentyp ist sehr einfach 
durchzufuhren, da der erste Zeitreihentyp klassifiziert wird, wenn der Graph des Informationsflusses eine in 
etwa gekrummte Form aufweist 

Weiterhin ist es vorteilhaft, das Verfahren fur eine Zeitreihe einzusetzen, die durch ein gemessenes Elektro- 
kardiogramm-Signal (EKG) zur Verfugung gestellt wird. Mit der Ermittlung stochastischer Korrelationen 
45 zwischen den Abtastwerten der Zeitreihe wird eine Klassifikation der Zeitreihe in ein Elektrokardiogramm-Si- 
gnal (EKG), welches ein Herz beschreibt gefahrdet ist bezQglich des pldtzlichen Herztodes sowie in ein Elektro- 
kardiograram-Signal (EKG) eines ungefahrdeten Herzens, moglich. Dadurch ist es mdglich, friihzeitig eine 
Gef ahrdung zu erkennen und eine Behandlung gegen den pldtzlichen Herztod einzuleiten. 

Anhand der Fig, 1 bis 5, die ein Ausfuhrungsbeispiel darstellen, wird die Erfindung im weiteren naher erlautert 
50 Es zeigen 

Fig. 1 ein Ablaufdiagramm, in dem das erfindungsgeraaBe Verfahren dargestellt ist; 

Fig. 2 ein Ablaufdiagramm, in dem die Weiterbildung des erdindungsgemaBen Verfahrens gemaB Patentan- 
spruch 3 beschrieben ist; 

Fig. 3 ein Blockdiagramm, in dem verschiedene mdgliche Zeitreihen, die klassifiziert werden, dargestellt sind; 
55 Fig. 4 ein Blockdiagramm, welches einen Rechner, der zur Durchf uhrung des erfindungsgemaBen Verfahrens 
notwendigerweise verwendet wird, dargestellt ist; 

Fig. 5 ein Diagramm, in dem qualitativ ein Graph eines ermittelten Informationsflusses fQr zukQnftige Werte 
fQr eine chaotische Zeitreihe, eine Zeitreihe, die nichtlineare Korrelationen zwischen ihren Abtastwerten auf- 
weist, sowie eine Zeitreihe, deren Abtastwerte stochastisch voneinander unabhangig sind, dargestellt sind 
60 In Fig. 1 ist dargestellt, daB in einem ersten Schritt des erfindungsgemaBen Verfahrens die Zeitreihe, welche 
eine vorgebbare Anzahl von Abtastwerten aufweist, gemessen wird 101. Die Messung erfolgt durch ein MeBge- 
rat MG, welches sowohl analoge oder digitale Signale miBt und einem Rechner R zufOhrt (vgL Fig. 4). Von dem 
Rechner R werden bedingte Entropien H(n|n— 1 ... 1) fQr die einzelnen Abtastwerte der Zeitreihe ermitteit 102. 
Verschiedene Vorgehensweisen zur Ermittlung der bedingten Entropien H(n|n— I ... 1) sind bekannt [2J [3} Im 
Rahmen dieses Dokumentes wird fflr die bedingten Entropien H(n|n - 1 ... 1) beispielsweise folgende Definition 
verwendet, welche jedoch nicht die Moglichkeit der Verwendung anderer Definitionen ira Rahmen des erfin- 
dungsgemaBen Verfahrens einschrankt: 
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HOJn - L . .1) = ^ ^ i) . log(p(j|i)) (1) 

5 

wobei mh 

H{n|n— 1 . . . 1) jeweils die bedingten Entropien bezeichnet werden, 

n eine Lange einer Sequenz beriicksichtigter Abtastwerte der Zeitreihe bezeichnet win!, 

kn (kn - m") eine Anzahl verschiedener Sequenzen beriicksichtigter Abtastwerte der Ltnge n bezeichnet wird, 10 
m eine Anzahl von Werten bezeichnet wird, die die Abtastwerte annehmen konnen, 
p(j, i) Verbundwahrscheinlichkeiten bezeichnet werden, und 
pG|i) bedingte Wahrscheinlichkeiten bezeichnet werden. 

Es ist bei dem erfindungsgeraaBen Verfahren vorgesehen, die bedingten Entropien H(n|n— 1 ... 1) far die 
vorgebbare Anzahl von Abtastwerten, die die Zeitreihe aufweist, zu bestimraen. Es ist jedoch ebenso vorgese- 15 
hen, einige bedingte Entropien H(n|n-1 ... 1) nicht zu ermitteln, und somit die entsprechenden Abtastwerte 
nicht zu berflcksichtigen. Dies entspricht einer Verringerung der Anzahl von Abtastwerten. Die Anzahl beriick- 
sichtigter Abtastwerte der Zeitreihe spiegeln direkt die Genauigkeit des erfindungsgemaBen Verfahrens bezGg- 
iich der Klassifikation der Zeitreihe wider. 

Die Anzahl von Werten m, die die Abtastwerte annehmen kdnnen, ist vorgebbar. Die Werte kdnnen, mussen 20 
jedoch nicht iiber konstante IntervaJIe verteilt sein. 

Es kdnnen ebenso fQr verschiedene Klassifikationen unterschiedliche mogliche Werte von Abtastwerten 
vorgegeben werden. Ein Satz von vorgebbaren Werten der Anzahl in wird im folgenden als eine Partition p 
bezeichnet Somit bezeichnet die Partition p einen Satz disjunkter IntervaJIe Bi d k 

25 

m 

P - ( B i)i =1 / (J B i = A ™d . B ± n Bj = 0 for i * j 

i=l >a^* 



30 



Hierbei bezeichnen i und j einen ersten Laufindex und einen zweiten Laufindex. 
Somit ergibt sich als eine Blockentropie 



k 35 

- EP 1 ' P (n) ■ logfpi'Pcn)) f5>S 
i = l V J 



Hierbei bezeichnet p"*(n) die Wahrscheinlichkeit des Auftretens eines Abtastwerts, der fur die Partition B den *° 
Abtastwert i aufweist bei einer Sequenz der Lange n. 
Eine Entropie f ur eine vorgebbare Anzahl zukunftiger Abtastzeitpunkte p ist gegeben durch 



H P (n,p) = - £ ZP i ' j ' P (n,p)-logfp i '^P (n , p) > ) (4) 



i = lj=l 



50 



Hierbei bezeichnet p y<p (n,p) die Verbundwahrscheinlichkeit des Auftretens eines Abtastwertes i fur die 
Sequenz der Lange n und das Auftreten des Abtastwertes j zu einem Zeitpunkt, der die vorgebbare Anzahl 
zukunftiger Abtastzeitpunkte voraus ist im Rahmen der Partition p. Eine bedingte Entropie, jeweils unter der 
Voraussetzung, daB der zeitlich direkt vorangegangene Abtastzeitpunkt bekannt ist, wird mit HH(n+ l)]n » 
bezeichnet r 55 

Ein InformationsfluB fur die vorgebbare Anzahl zukunftiger Abtastzeitpunkte p fQr eine bestimrate Partition B 
wird gebildet nach folgender Vorschrift: 

l£ = lim lP(n + p, n + ljn— 1) (6) 60 

n— »oo 

Hierbei ergibt sich l\n + p,n + 1 |n . . . 1) aus: 
I*(n + pji + 1 |n . . . 1) - H**(n + p|n . . . 1) + H**(n + p|n + 1 . . . 1) (7) 65 
Die Partition p ist als eine infitesimaJe Partitionierung definiert, so daB gilt: £ « diameter (P) 0, wobei mit 
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diameter (B) eine jeweils grdBte Zellinge bezeichnet wird. m . r , 

Der IiuSrmationsfiuB Pp ist aUgemein ein MaB fur die statistisch nichtlineare Korrelati n zwischen zwei 
freien Variable* in diesem Verf ahren zwischen den Abtastwerten verschiedener Abtastzeitpunkte der Zeitreihe. 

Der InformationsfluB tf» p wird somit bei dem erfindungsgemiBen Verfahren fur eine vorgebbare Anzahl 
zukunftiger Abtastzeitpunkte p abhingig von einer vorgebbaren Anzahl vergangener Abtastwerte n, die die 

^Au! Sen bedhSen EntSpien wird in einem dritten Schritt 103 mindestens ein InformationsfluB tfp bestimmt- 
Eto Graph^Funktion des Informationsflusses Pp weist fur unterschiedUche Zeitreihen unterschiedhche 

C ^„^tiS^ der InformationsfluB I* P einer Partition B in einer idealen Niherung 

einen konstanten. waagrechten Verlauf Ober den Abtastwerten p auf. . ... _ „ „io*;~- 

Quahtativ ergibt sicfa fflr den InformatinsfluB tf> P einer Zeitreihe, deren Abt^twerte nichthneare Korrelatio- 
nen aufweisen. eine monoton faUende, parabelahnlich gekrQmmte Funktion ZT1. Dies entspncht einem ersten 

Z wSS l jedcTdie Abtastwerte untereinander keinerlei Korrelationen auf, so ist qualitativ ein stett, in etwa 
linear abfallender Graph des Informationsflusses I» p for zukflnftige Abtastwerte gegeben. Dies ist ersichthch 
S die l^eriegung, daB bei nicht vorhandener Korrelation zukQnftiger Abtastwerte in keinster Weise 
vorhergesagt werden kdnnen. und somit auch keinerlei Information Ober zukflnftige Abtastwerte vorhanden 
sind Dies ist eben fur eine Zeitreihe, deren Abtastwerte nichtlineare Korrelationen aufweisen, nicht derFalL 
In einem letzten Schritt 104 wird anhand des Infonnationsflusses F P erne Klassifikation durehgefOhrt. Diese 

laassifikationkannjenachAnwendungsbereichunterschiedlicherArtsein. 

Eine sehr einfach Klassifikation, die sich jedoch fflr einige Arten von Zeitreihen als eine vorteilhafte und 
ausreichendeWeiterbadungdesVerfahrensemeist,Uegtinemer^maren"KlassuTkaOon. 

Hierbei wird anhand des Graphen des Informationsflusses Pp for zukunftige Abtestwerte 201 in einem 
Prflfungsschritt 202 Oberpruftob der Graph in etwa gekriiiiimtist,oderobersteaimearfallt(vgLFig : 2^ 

Weist die Form des Graphen eine parabeiahnUche, leicht gekrflmmte, faUende Form auf, so wird die Zeitreihe 
als der erste Zeitreihentyp ZT1 klassifiziert. Dies entspricht bei einer Zeitreihe, die durch e in ' gemessenes 
Kardiogramm-Signal (EKG) gegeben ist, einer Klassifikation des Elektrokardiogramm-Signak (EKG) in em 
Eektrokardiogramm-Signal (EKG) eines gef ihrdeten Herzens bezflglich des Pldtzjichen Herztodes. 

Weist jedoch der Graph eine steU abf allende. lineare Form auf, so wird die Zeitreihe in den zweiten Zeitreihen- 
typ ZT2 klassifiziert 203. Dies entspricht fflr das Beispiel des Hektrokardiogramm-Signals der Klassifikation des 
EKG-Sienals als ein EKG-Signal eines ungef ahrdeten Herzens bezugUch des pldtzlichen Herztodes. 

In Fig! 3 sind verschiedene M8gUchkeiten fur die Arten einer Zeitreihe, fur die das Verfahren einsetzbar ist, 
angegeben 301. Diese Aufzahlung weist jedoch keinerlei einschrankenden Charakter auf. Das Verfahren ist fflr 
iede Art von Zeitreihe verwendbar, bei dem es gilt, nichtUneare Korrelationen zwischen den Abtastwerten der 
Zeitreihe zu ermitteln, und die Zeitreihe anhand dieser nichtlinearen Korrelationen, die sich ui dem Informa- 
tionsfluB widerspiegeln, zu klassifizieren. 
Die Zeitreihe kann beispielsweise sein: 

— ein Elektrokardiogramm-Signal (EKG) 302, 

— einElektroencephalogramm-Signal(EEG)303, .... 

— ein Signal welches den Verlauf der Sauerstoff spannung eines Gehirns beschreibt 304. 

In Fig. 4 ist der Rechner R dargesteUt, mit dem das erfindungsgemlBe Verfahren notwendigerweise durchge- 
45 f0 ^^jJ :|iiier R verarbeitet dig von dem MeBgerat MG aufgenommenen, und dem Rechner R zugefflhrten 

ZC Hierbe?ist es nicht von Bedeutung, ob die BUdung der Abtastwerte aus dem m6glicherweiseanalogen Signal 
in dem MeBgerat MG oder in dem Rechner R durehgefOhrt wird. Beide Varianten sind fur das erfindungsgemaBe 
so Verfahren vorgesehen. _ v . _ . . , , .-——v 

Das MeBgerat MG kann beispielsweise ein Elektrokardiograph (EKG), ein Elektroencepahlograph (EEG) 
oder auch ein GerSt sein, welches nach dem in [2] dargesteUten Verfahren arbeitet 

Das Klassifikationsergebnis, welches durch den Rechner R auf die im vongen beschnebene Weise erauttelt 
wird, wird in einem Mittel zur Weiterverarbeitung WV weiterverarbeitet, beispielsweise einem Bemitzer darge- 
steUt Dieses Mittel WV kann beispielsweise ein Drucker, ein BUdschirm oder auch ein Lautsprecher sein, uber 
das ein akustisches oder visueUes Signal an einen Benutzer weitergegeben wird. 

ImRahmen dieses Dokuments wurden folgende Verof f entUchungen zitiert: ien!oa . 
r 11 G. Morfill, Komplexitatsanalyse in der Kardiologie, Physikalische Blatter, VoL 50. Nr. 2. S. 156 bis 160 1994; 
[2] W. Ebeling et al Entropy, Transinformation ans Word Distribution of Information-Carrymg Seqences, 
International Journal of Bifrucation and Chaos, VoL 5, Nr. 1,S. 51 -61, 1995; 

[3] D. Wolpert et al, Estimation Functi ns of ProbaUty Distributions from a finite Set of Samples, Physical 
ReviewE,VoL52,Nr.6,S.6841-6854,Dezemberl995; . . 

[4] UCOX. GMS, GeseUschaft fflr Medizinische Sondentechnik mbH, Advanced Tissue Momtonng; 
[5]DE39 12 028A1. 

PatentansprQche 

1. Verfahren zur Klassifikati n einer meBbaren Zeitreihe, die eine vorgebbare Anzahl von Abtastwerten 
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aufweist, beispielsweise eines elektrischen Signals, durch einen Rechner, 

— bei dem bedingte Entropien ermittelt werden, 

— bei den aus den bedingten Entropien mindestens ein InfonnationsfluB for eine vorgebbare Anzahl 
zukOnftiger Abtastzehpunkte (p) bestimmt wird, und 

— bei den anhand des Informationsflusses eine Klassifikation der Zeitreihe durchgefOhrt wird 

2. Verfahren nach Anspruch 1, bei dem alle bedingten Entropien der Abtastwerte der Zeitreihe ennittelt 
werden. 

3. Verfahren nach Anspruch 1 oder 2, bei dem die bedingten Entropien ermittelt werden nach der Vorschrift 

k n -l m 

H(njn — L . .1) = — Z Z p( > i) - log(p(j]i)) 
i=lj=l 

wobei mit 15 

H(n|n - 1 . . . 1) jeweils die bedingten Entropien bezeichnet werden, 

n eine Lange einer Sequenz berflcksichtigter Abtastwerte der Zeitreihe bezeichnet wird, 

ko (kn - m n ) eine Anzahl verschiedener Sequenzen berOcksichtigter Abtastwerte der Lange n bezeichnet 

20 

m eine Anzahl von Werten bezeichnet wird, die die Abtastwerte annehmen kdnnen, 
p(j, i) Verbundwahrscheinlichkeiten bezeichnet werden, und 
pOI 0 bedingte Wahrscheinlichkeiten bezeichnet werdea 

4. Verfahren nach einem der AnsprOche 1 bis 3, bei dem der InfonnationsfluB (frn + 1 ; n + pin ... 1)) fur eine 
vorgebbare Anzahl zukflnftigerAbtastzeitpunkte(p) bestimmt wird nach der Vorschrift 25 

I^n + l;n + p|n... 1) - H^n + p|n... 1) - H^n + p|n + 1...1) 

wobei 

P eine Partition eindeutig bezeichnet durch die eine Anzahl m von Werten, die die Abtastwerte annehmen 30 
konnen,vorgegebenwird 

5. Verfahren nach einem der AnsprOche 1 bis 4, bei dem bei der Klassifikation die Zeitreihe entweder in 
einen ersten Zeitreihentyp, wenn der Graph des Informationsflusses eine gekrOmmte Form aufweist, oder in 
einen zweiten Zeitreihentyp klassifiziert wird, wenn eine gekrOmmte Form des Graphen nicht vorliegt. 

6. Verwendung des Verfahrens nach einem der AnsprOche 1 bis 5, wobei als Zeitreihe ein Elektrokardio- 35 
gramm-Signal (EKG) gemessen wird. 

^Verwendunjg : des i Verfahrens nach einen der AnsprOche 1 bis 5, wobei als Zeitreihe ein Elektroencephalo- 
gramm-Signal (EEG) gemessen wird. 

8 ^ e T e ^ dUr ? dM Verfahr f m , n ? ch einem der AnsprOche 1 bis 5, wobei als Zeitreihe ein Signal gemessen 
wird, das den Spannungsverlauf ernes Gehirndrucks beschreibt 
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